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Abstract 

Using a relativistic nuclear optical potential consisting of a Lorentz scalar, Vg, and the time-like 
component of a four- vector potential, Vq, we calculate elastic scattering differential cross sections 
and polarizations for p+'^He at intermediate energies for which experimental data are available. 
We also calculate the differential cross sections and analyzing powers for p-|-^'^He at intermediate 
energies and compare with the few available experimental data. 

PACS numbers: 24.10.-i,24.10.Eq,24.10.Jv 
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I. INTRODUCTION 



The scattering of medium-energy nucleons from nuclei can provide information about 
both nuclear structure and the A^A^-interaction. In particular, the study of proton-nucleus 
scattering at intermediate energies is a useful method for determining accurate nuclear 
matter distributions in stable nuclei [|T|. Recently, it has also been used extensively to 
determine the extension of the nuclear matter density of exotic nuclei. In particular, in 
a recent set of experiments, inverse kinematics was used to study the elastic scattering of 
p_l_4,6,8jjg determine their matter densities |]^, |^. 

Elastic scattering at intermediate energies can be calculated theoretically using either 
a non-relativistic @, |^ or a relativistic optical [H, model, usually with about the same 
good results. An advantage of a relativistic optical model is that its two potentials si- 
multaneously determine both the central and spin-orbit interactions. When the relativistic 
impulse approximation (RIA) is used [§, these interactions may be determined in terms of 
the corresponding scalar and vector nuclear densities @. 

In this work we use a relativistic nuclear optical potential constructed from a Lorentz 
scalar, Vs, and the time-like component of the four-vector potential, Vq, to calculate elastic 
scattering angular distributions and analyzing powers for p+^He at intermediate energies 
and compare two different fits to the p+'^He data and also to the RIA results. We also 
calculate angular distributions and analyzing powers of the elastic scattering of p+^'^He and 
compare an extension of the more physical of our two fits with the RIA and the few available 
experimental data. 

II. THE DIRAC OPTICAL POTENTIAL 

In relativistic optical model analyses of intermediate energy scattering the Dirac 
equation is usually used in the form, 

{a ■ p + (3[m + Vs{r)] + [Vo{r) + K(r)]}^(r) = E^{f), (1) 

where Vg is an attractive Lorentz scalar potential, Vq is the repulsive time-like component 
of a four-vector potential and is the Coulomb potential. The choice of the potentials is 
motivated by meson exchange considerations and simplicity. The simplest meson exchange 
interaction possessing a certain justification on physical grounds, which is also capable of 

2 



providing nuclear saturation properties, takes into account the exchange of an attractive 
intermediate range isoscalar scalar meson and a repulsive shorter range isoscalar vector 



meson |T0|, |Tl|, |T^. The two corresponding potentials K and Vq play an essential role 
in the description of both elastic scattering and polarization data at intermediate energies 
since their sum is the principal contribution to the central potential while their difference 



determines the spin-orbit interaction |]T3[. The scalar, Vg, and vector, Vq, optical potentials 



used in the analyses here can be written as 

Vs = UJuAr)+^WsfwAr) 

Vo = Uofuo{r) + iWofwoir) (2) 

with each of the potentials possessing both real and imaginary parts with possibly different 
radial dependences. 

Based on nuclear structure calculations, we take the radial dependence /(r), for the case 
of p+'^He, to have a Gaussian form 

/(r) = exp[-rVro]. (3) 

III. RESULTS 

A. p+^He elastic scattering 

Relativistic Hartree calculations of the ^He nucleus yield almost identical scalar and vector 



densities but an rms radius about 30% larger than the experimental one [Q. To obtain a 
physically reasonable RIA potential for this system, we use equal scalar and vector matter 
densities having the rms radius of 1.49 fm found in Ref. M together with the relativistic 



Love-Franey NN t-matrix of Ref. |TB|. The results, represented by full lines in Figs. (|Tp-(^, 
show that the RIA agrees well with the experimental data at low momentum transfer but 
deviates substantially from the experimental results as the momentum transfer increases. 
We also note that the RIA angular distributions do not reproduce the oscilatory structure 
seen in the experimental data, decreasing monotonically instead. The polarizations are 
not reproduced either, except at very low values of the momentum transfer. Due to these 
differences, we have tried to fit the experimental data of Ref. using a Dirac optical 
potential of the form given in equations 



In a first attempt, labeled with the dashed lines in the figures, the parameters were left 
free to vary so as to obtain the simultaneous best fit to the experimental angular distributions 
and polarizations at the three values of the laboratory energy, Eiab- The same geometrical 
parameters were used at the three energies, while the strengths were assumed to vary linearly 
with the laboratory energy as V = Vq + Vi* Eiab, a common parametrization of optical model 
strengths. The best fit parameters are given in Table |. As can be seen in the figures, we 
indeed obtain a fairly reasonable fit to the angular distribution and polarization at all three 
energies. However, an interpretation of the imaginary part of the potential in the context 
of a RIA potential would require a A^A^ total cross section about 2.5 times the physical 
one, which makes the fit unsatisfactory on physical grounds. The radii of the imaginary 
potentials are in agreement with the matter radius of Ref. 0], while the radii of the real 
potentials are found to be about 10% bigger. These also disagree with what one would 
expect from a RIA potential, which usually yields real radii close to the matter radius and 
imaginary ones slightly larger. 

To obtain a more physically reasonable fit to the p+^He data, we reduced the imaginary 
strengths of the first fit by a factor of 3, reset the radii to the value of tq = 1.22 fm, 
corresponding to an rms radius of 1.49 fm 0, and then let the parameters vary freely 
once again. The results, labeled with dotted lines in the figures, present slightly poorer 
agreement with the experimental data, but continue to reproduce the oscillations in the 
angular distributions that are not obtained in the RIA calculation (full lines). The imaginary 
potentials of the second fit are consistent with those of the RIA potential, corresponding to 
a physically reasonable A^A^ total cross section. The radii are also closer to those of a RIA 
potential. 
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for fits to the p+^He experimental data. 
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FIG. 1: Differential cross section of p+^He, as a function of the momentum transfer, calculated 
for Eiah = 0.561, 0.800 and 1.029 GeV. The solid lines represents the RIA calculation while the 
dashed and dotted ones stand for the first and second fits. The experimental data are labeled with 
diamonds. 

B. p+^'^He elastic scattering 

We next considered the elastic scattering of p+^'^He at laboratory energies of Ei^b = 0.721 
GeV and Eiab — 0.678 GeV, respectively. To perform the RIA calculation for ^'^He, we 
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FIG. 2: Polarization of p+'^He, function of the momentum transfer, calculated for Eiab = 0.561, 
0.800, 1.029 GeV. The solid lines represents the RIA approximation while the dashed and dotted 
ones stand for the first and second fits. The experimental data are labeled with diamonds. 

followed the same procedure used for ''He. We generated equal scalar and vector matter 
densities using the appropriate Gaussian harmonic oscilator (HO) geometrical parameters 
of Ref . , together with the relativistic Love-Franey NN t-matrix of Ref . . The results 
of the RIA calculation in Figs. (0) and agree well with the experimental data in the 
range over which the data exist. We emphazise that the RIA also describes the p+^He 
data well in this range of momentum transfers but, at larger momentum transfers, deviates 
substantially from the data in form and absolute values. 

For the cases of p+^'^He, the data is insufficient to attempt a fit. We have thus attempted 
to extend our fits to the ^He data by interpreting them in the context of a simple RIA 
calculation. To do so, we extract the effective strengths, Vq, = VFj)[7rrg J^/^, that, 

when multiplied by the ^He density, would yield the potentials obtained in our fits to the 
j9+^He data. We have multiplied these strengths by the appropriate ^'^He densities, obtained 
using the Gaussian HO parametrization of Ref. P] to obtain effective Dirac potentials for 
these systems. 
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where ( = 2,4 for ^He and ^He, respectively. The radii were modified shghtly to take into 
account the differences between the experimental ^He radius of Ref. and those obtained 
in the fits. 

Our calculations for these two nuclei do not reproduce the experimental data as well as 
the RIA calculation does, as can be seen in Figs. @ and (^), which leads us to conclude 
that the parameters obtained in the p+^He fit cannot be so easily extended to the p+^'^He 
elastic scattering. However, we feel that they can transmit an idea of the discrepancies with 
an RIA calculation that the experimental values could show. In Figure (^) we show the 
polarizations expected for proton scattering from these nuclei, for both the RIA and our 
extension of the ^He fits. The differences between the two are again quite large. However, 
no data exist in this case. By comparision with the results for ^He it is difficult to say 
what one could expect of the experimental angular distribution and polarization in p+^'^He 
scattering at larger values of the momentum transfer. 




FIG. 3: Differential cross sections of p+^'^He at E^^ab = 0.721, 0.678 GeV, respectively. The dashed 
hne corresponds to the RIA approximation while the dotted one represents the second fit. The 
experimental data are labeled with diamonds. 
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FIG. 4: Polarization of p+^'®He at Eiab = 0.721, 0.678 GeV respectively. The dashed line corre- 
sponds to the RIA approximation while the dotted one represents the second fit. 

IV. CONCLUSION 

We have calculated p+^'^'^He elastic scattering differential cross sections and polarizations 
using the relativistic impulse approximation and an adjusted Dirac optical potential. We 
have shown that the RIA results obtained using the parameters of Ref. describe all three 
systems well at low values of the momentum transfer. However, in the case of p+^He, for 
which more data exist, the RIA deviates substantially from the data at higher momentum 
transfers. As an alternative to the RIA we have ajusted the parameters of a Dirac optical 
potential and then attempted to extend the parameters obtained to the cases of p+^'^He, 
but with unsatisfactory results. Based on our analyses it is difficult to say what one could 
expect at higher values of transferred momenta in the cases of p+^'^He. More data are 
needed to elucidate this situation. 

In closing, we should mention that a possible reason for the failure of RIA in describing 
the p+^He data is that it has been pushed beyond the limit of its applicability in the simple 
manner in which it is used here. The decomposition of the effective Dirac potential into a 
scalar potential and the fourth component of vector is valid in the target rest frame. The 
scattering calculation, however, is performed in the CM frame and the boost that carries 
one frame to the other will convert the vector fourth component potential into a full vector 
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potential. For proton scattering on a system such as Ca the vector components introduced 
are small and can be neglected. For extremely light systems, such as those studied here, 
this is not the case. We plan to examine the importance of this effect in the future. 
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